Weak Localization in Metallic Granular Media 
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We investigate the interference correction to the conductivity of a medium consisting of metallic 
grains connected by tunnel junctions. Tunneling conductance between the grains, e^g^j-nh, is 
assumed to be large, ijt > 1- We demonstrate that the weak localization correction to conductivity 
exhibits a crossover at temperature T ~ g%&, where 5 is the mean level spacing in a single grain. At 
the crossover, the phase relaxation time determined by the electron-electron interaction becomes of 
the order of the dwell time of an electron in a grain. Below the crossover temperature, the granular 
array behaves as a continuous medium, while above the crossover the weak localization effect is 
largely a single-junction phenomenon. We elucidate the signatures of the granular structure in the 
temperature and magnetic field dependence of the weak localization correction. 

PACS numbers: 73.23.-b, 72.15.Rn, 73.23.Hk 



I. INTRODUCTION 

Quantum effects in conduction of two-dimensional dis- 
ordered electron systems draw attention of both experi- 
mentalists and theorists for decades. The interest is moti- 
vated in part by the interplay between several fundamen- 
tal physical phenomena, such as quantum interference, 
localization, superconductivity, and single-electron tun- 
neling occurring in these systems. The interplay affects 
the properties of normal conductors' 2 ^, nominally su- 
perconducting films 4 ^, and arrays of junctions^. Quan- 
tum effects become increasingly important at sheet con- 
ductances decreasing towards the fundamental quantum 
value of Gq = e 2 /irh. The interpretation of some of the 
most intriguing data, however, may depend on whether 
the investigated conductors are homogeneous or granu- 
lar. While this question has a definite answer in the case 
of an array£ of lythographically defined junctions, it is 
less clear for nominally homogeneous deposited metallic 
films 4 ^ or electron gases in semiconductor heterostruc- 
tures^. Checking the samples homogeneity traditionally 
involves application of auxiliary techniques, such as local 
probe spectroscopjs^i. 

We demonstrate that information about the granular- 
ity of a conductor is contained in the temperature and 
magnetic field dependence of the weak localization (WL) 
correction to the conductivity. The granular structure of 
a conductor affects the correction even at high film con- 
ductivity, cr ^> Gq. While being universal at the lowest 
temperatures and magnetic fields, the WL correction be- 
comes structure-dependent at higher values of field and 
temperature. The corresponding crossover temperature 
is of the order of (cto/Gq) 2 ^ where the mean-level spac- 
ing 6 in a single grain is inversely proportional to the 
grain volume. The field dependence of the WL correc- 
tion at low temperatures exhibits two crossovers. These 
are associated with a significant change in structure of 



closed electron trajectories, allowing for phase-coherent 
electron motion. 

The WL correction in a homogeneous medium origi- 
nates from the quantum interference of electrons moving 
along self-intersecting trajectories^ and is proportional to 
the return probability of an electron diffusing in a disor- 
dered medium. In one- or two-dimensional conductors 
this probability diverges due to the contribution com- 
ing from long trajectories. For a fully coherent electron 
propagation, this divergence would lead to a divergent 
WL correction. Finite phase relaxation time makes 
sufficiently long trajectories ineffective for the interfer- 
ence, and limits the correction. In a two-dimensional 
conductor, the WL correction to conductivity is 8a = 
— (Gq/2t:)Iti(t^/t), where r is the electron momentum 
relaxation time. There are various mechanisms of the 
electron phase relaxation, some of them being material- 
specific 8 . The most generic and common for all the con- 
ductors mechanism stems from the electron-electron in- 
teraction^. It yields l/r^ ~ T^Gq/ctqJI^Gq/o-o) and 
provides the temperature dependence of the conductiv- 
ity, 



a Q -(G Q /2n) ln(T*/T) 



(1) 



with T* ~ (Jq/Gqt. The typical area under an electron 
trajectory which barely preserves coherence, L\ — Dt^, 
depends on the electron diffusion constant. Magnetic 
field B significantly affects the WL correction if the cor- 
responding magnetic flux through a contour of area L 2 , 
exceeds the quantum $o- This makes the magnetoresis- 
tance measurement a useful tool for the investigation of 
the electron interference. 

To model a granular medium, we consider a regular 
two-dimensional array of grains of size d connected by 
tunnel junctions. The grains have internal disorder, but 
are characterized by conductance far exceeding the con- 
ductance 3tGq of a single tunnel junction. The classi- 
cal conductivity of a square array is thus ctq — GqQt- 



2 



It corresponds^ to the effective electron diffusion con- 
stant D = 7r^ 1 gr r Sd 2 . In the absence of phase relax- 
ation, an electron may pass through any number of junc- 
tions coherently. It will result in a diverging WL cor- 
rection, just like in a homogeneous conductor. The 
electron-electron interaction limits the phase relaxation 
time, yielding 1/r^ ~ T/g^. As long as the correspond- 
ing length Ltf, ~ dw g^,8/T significantly exceeds d, an 
electron may return to a grain coherently after passing 
many junctions, and the inhomogeneity of the granular 
medium is irrelevant. The comparison of with d de- 
fines a crossover temperature, 



(2) 



Roughly, above the crossover temperature the electron 
trajectories contributing to the WL do not cross more 
than a single junction. In this regime granular medium 
behaves similarly to a single grain connected to highly 
conducting leads by tunnel junctions of conductance gr- 
The WL correction at T > T CI comes from electron 
trajectories that pass through a single tunnel junction. 
Electrons moving along longer trajectories, which include 
more junctions, have much smaller probability of a phase- 
coherent return. We find that already the shortest inter- 
grain trajectories (see Fig. 2 in Section Hvjl provide the 
temperature dependence of the WL correction, 



WL 



-A- 



(3) 



with A being a geometry-dependent coefficient of order 
one. In deriving Eq. J3J|, we assume that <7t is much 
smaller than the number of channels in the inter-grain 
tunnel junction, although qt 1. 

Equation does not account for the phase relaxation 
rate within the grains. At a sufficiently high tempera- 
ture T > T* the latter exceeds the electron escape rate 
gi>5 from a grain, which leads to a suppression of the 
WL correction below the value Eq. |J3}. The charac- 
teristic scale T* here depends on the intra-grain phase 
relaxation mechanism. Assuming that it is due to the 
electron-electron interactional, and that the dimension- 
less conductance of the grain g gr is large, <? gr > g^ 7 we 
find 



9 & nr 



(4) 



In a more exotic case of a smaller grain conductance, 
9t ^ 5gr ^ ffT, the temperature T* still exceeds sig- 
nificantly T cr , but the specific relation between the two 
temperature scales depends on the grain shape, and is 
different for disk-like or dome-like grains. 

We turn now to the discussion of the magnetic field 
effect on the weak localization in the granular medium. 
To determine the characteristic field suppressing the in- 
terference correction to conductivity, we need to estimate 
the directed area covered by a typical closed electron tra- 
jectory^. For a single grain, such area is of the order 








FIG. 1: A typical diffusive trajectory in a granular array. The 
directed area Seg consists of two components. The first one 
is the combined contribution of separate grains, see the first 
term in Eq. @. The second component is the area under 
the coarse-grained trajectory, which is the counterpart of the 
directed area under an electron trajectory in a homogeneous 
disordered sample. 



d 2 \fg~&Jg~T and is limited by the electron dwelling time. 
At low temperatures, T <C T cr , the number of grains vis- 
ited by an electron having a typical closed trajectory, is 
of the order of T cr /T. The single-grain directed areas 
have random signs, so the estimate for the full directed 
area is 



S c fi 



' -*cr<7gr 

Tgr 



(•5) 



The first term here corresponds to the sum of the directed 
areas under the electron trajectory within the grains vis- 
ited by electron; the second, conventional^, term comes 
from the fact that the "visited" grains form a closed con- 
tour of an area IA, see Fig. 1. The characteristic level 
of the field necessary to affect (5<7wl is found from the 
condition S c ftB<p ~ $o- We see now, that even within 
the temperature range T < T cr , the granularity of the 
material matters. 

At the lowest temperatures the characteristic field co- 
incides with that of a film with the corresponding value 
of diffusion coefficient, 

B *~%f-> T<<Tj f- (6) 

Uj j. cr yg r 

At higher temperatures, the characteristic field is 



T CI — «T« T cr . (7) 

3gr 



The higher the applied field, the shorter are the trajecto- 
ries contributing to the interference correction, and the 
smaller the correction is. Such trajectories span only a 




3 




single grain provided the field B is of the order or higher 
than 



(8) 



At B > £>^ s , even the single-grain Cooperon is sup- 
pressed. Consequently, Eq. JSJ defines the characteristic 
field suppressing at T > T cr the WL correction Eq. J3J), 
which stems from the transitions within the closest grain 
pairs. 

To develop a quantitative theory of the interference 
correction, we derive the expression for the weak local- 
ization correction and adapt the Cooperon equation for 
granular medium. 



II. CONDUCTANCE AND WEAK 
LOCALIZATION IN A GRANULAR ARRAY 

Tunnel junctions between metallic grains are described 
adequately by a model with an infinitely large number of 
weakly-conducting channels. Within this model, one can 
use the tunneling Hamiltonian formalism for evaluation 
of the conductivity of the granular array. In this formal- 
ism, tunneling between the grains m and n is described 
by the Hamiltonian 



after 



tape 



iV mn trf 



ai a a Pa + h.c. =y~] / dri (9) 

„ J m 



x / dr 2 t(r 1 ,r 2 y v ^ t ^ m(J (r 1 )i, na {r 2 ) + h.c. 



where the points r\ and r 2 belong to the grains m and 
n, respectively, V mn is the voltage applied to the barrier, 
a € m and (3 £ n are exact single-particle states, and a 
is the spin index. In the limit of thin barrier, the tunnel 
amplitude t significantly deviates from zero only if the 
vectors r\ and r 2 refer to two closest to each other points 
at opposite sides of the interface, 



t(r 1 ,r 2 )=a6(y-y')d z d z ,6(z)5{z'). 



(10) 



Here the coordinate y runs along the interface S, and 
transverse coordinates z in the grain m and z' in the 
grain n are defined in such a way that at the interface 
z = z' = 0. (We wrote Eq. I|10f> for the planar geometry, 
generalization to three-dimensional arrays is straightfor- 
ward). The constant a is of the order of magnitude of 
VT /vkF, where v is the electron density of states of the 
material of the grains, and T is the transmission coef- 
ficient of the barrier. The numerical factor in a can be 
related to the measurable quantity, the barrier conduc- 
tance <7t- Using Eq. (|10fl . one may express the tunnel 
amplitude in terms of the eigenfunctions \a and xp °f 
an electron in the grains m and n, respectively (see, e.g., 
Ref.lH, 

t af} = a dyd z x a {y, z )dz<Xf}(y,z')\ z=z >=Q . (11) 
Js 



The current through the contact is defined as I = 

—eN m = —ie[HT, N m ], where N m is the number of par- 
ticles in the grain m, 

N m = 2J o) aa d aa . 

aa 

Calculating the average current through the barrier, we 
obtain 

J(i) = -e Re J dr x dr 2 {t(r 1 .r 2 )G^ m (r 2 ,r 1 ,t,t)( 

-r(r 1 ,r 2 )G^„(r 2 ,r 1 ,i,t)e- iy - t } , 

where G K is the Keldysh component of the matrix 
Green's function, and the subscripts m and n are intro- 
duced for convenience, in order to indicate which grain 
points Ti and r 2 belong to. We need now to calcu- 
late the function G K by perturbation theory in tunnel- 
ing Hamiltonian. Let us first discuss the first order and 
calculate the average conductance. Using the standard 
technique^, we obtain for the current in the frequency 
representation (terms which do not depend on the time 
difference would correspond to the Josephson effect and 
thus are dropped) 

f dio 
J(w) = 2eRe / dr x ...dr A — t*(n,r 2 t(r 3 ,r 4 ) (12) 

x Tr T x G m (ri,r 3 ,uj + eV mn )Gn{r4,r 2 ,u!) 

where t x is the Pauli matrix in the Keldysh space, G r , 
G„ n , and we use the standard representation 



G 



G R G K 
G A 



In the linear regime, it suffices to use the equilibrium 
function here, G K (E) = tanh(E/2T){G R (E) - G A {E)). 
Expressing the Green's functions in terms of the exact 
eigenfunctions, calculating the energy integrals, and sub- 
stituting the transmission amplitudes I|1U|) . we obtain 
/ = uV mn for the inter-grain current, and a — (e 2 /ir)gT 
for the Drude conductivity of the granular array. The 
introduced here dimensionless inter-grain tunneling con- 
ductance is 



9T = 47r 2 |a| 2 



dydy' 



(13) 



= A\a\ 2 / dydy' [d z 8A™ (G R (yz;y'z')) \ z=z , =0 \ , 
Js 

where £ a are the exact energy eigenvalues measured from 
the Fermi level in a grain, and angular brackets mean 
impurity averaging within a grain (the eigenfunctions in 
different grains are not correlated). In the last equation, 
(G R ) is the impurity-averaged Green's function evaluated 
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m n 




FIG. 2: Second-order correction to the conductivity. Black 
circles represent the tunnel amplitudes, and dashed lines de- 
note impurity scattering inside the grains 



at the Fermi energy. It is represented as the density of 
states v multiplied with a dimensionless function rapidly 
decaying with the distance y — y' . The characteristic 
length of that decay is given by the Fermi wavelength, 
and the integral in Eq. i|13|) is converging rapidly. There- 
fore the dimensionless function oiy — y' may be evaluated 
within the free-electron approximation' 1 ''. The precise 
shape of this function is not important for our purposes. 
Equation thus relates the tunnel conductance to the 
previously introduced constant a. 

We proceed now with the evaluation of weak lo- 
calization correction. The next-order contribution to 
the current (Fig. 2) contains four tunnel amplitudes 
and four Green's functions with the Keldysh structure 
Tr ( r, GGGG), where G is the matrix Green's func- 
tion in the Keldysh space. The trace of a product 
of several Green's function can only have the follow- 
ing structure: several first functions are retarded, fol- 
lowed by one Keldysh function and then a number of 
advanced functions. Thus, we have the combination of 
the type G K G A G A G A + G R G K G A G A + G R G R G K G A + 
G R G R G R G K . However, the second and third terms in 
this combination are considerably greater than the other 
two, since the impurity scattering inside the grains is the 
most effective if the impurity line connect advance and 
retarded, advanced and Keldysh, and Keldysh and re- 
tarded Green's functions, but not two retarded or two 
advanced ones. Thus, retaining only these two term o 9 ' 15 , 
we express the weak localization correction in terms of 
the Cooperon C mn in the time representation, 

2e 2 f 

Sctwl = — —Re / drx . . .dr 4 t*(r 1 ,r 2 )t(r 3 ,r 4 ) 

dtC mn {r 1 ,r i ;r 3 ,r 2 ;t,-t) , (14) 

-oo 



where the subscripts of the Cooperon indicate that it 
starts and ends in the grains to and n, respectively. Note 
that due to the structure of the tunneling amplitudes 
t(r,r'), point r± is just across the barrier from point r 2 
and similarly point r 3 is across the barrier from point r 4 . 
The Cooperon C can be presented in the form 



C m n(ri,r i ;r 2 ,r 3 ) = -Im (G R (n - r 4 )) 
v 

xlm{G R (r 2 ~ r 3 )) C mn (ri,r 3 ) . 



Rapid decay of functions (G R ) with the distance between 
the corresponding arguments makes points in pairs r± , r 4 
and r 2 , r 3 in the spatial integral of Eq. (|14|) to be within 
the Fermi wavelength from each other. On the other 
hand, the Cooperon C(ri,r 3 ) is generally a long-range 
function. Provided we are interested in times long com- 
pared to the intra-grain diffusion time, C almost does 
not change while its coordinates vary within respective 
grains. However, C mn (ri,r 3 ) with m ^ n may dif- 
fer significanlty from the value of single-grain Cooperon 
(to = n). Substituting this coarse-grained Cooperon 
Cmn into Eq. ijTHl and taking into account Eqs. ijTU|) 
and (13) i we obtain 



5 a 



WL 



e 2 ffT 
2ttv 2 



Re 



dtC mn (t, —t) 



(15) 



with to and n being the neighboring grains. Note that 
Eq. (|15f) is valid for any dimension, not just in 2D. 

The form H15|l of weak localization correction is valid 
provided the phase coherence between the grains barely 
survives, and C mn -C C nn at to 7^ n. This limit is real- 
ized at a sufficiently high temperature, T ^> T cr . Note 
also that the performed derivation, unlike the consider- 
ation of, e.g., Ref. assumes the limit of large number 
of channels taken at fixed value of c/t- 

To consider phase relaxation in a granular array, we 
derive now the proper equation for C nm in a granular 
medium. 



III. COOPERON IN A GRANULAR ARRAY 



Cooperon describes the probability amplitude of elec- 
tron return and in the case of a homogeneous medium 
with electron diffusion coefficient D obeys the equation 



5 



at 



-D 



JL-ilA(r,-t/2)-i-A(r,t/2) 
or c c 



1 2 



C(r,r';t,t') = 5 (r - r') 5(t - t'). 



(16) 



Here the vector potential A accounts for the fluctuating 
electric fields representing the effect of electron-electron 
interactions, and should be considered as a Gaussian clas- 
sical random variable with zero average. 

In order to adapt Eq. (jTB)l to the case of a granular 
medium, it is convenient to perform a gauge transforma- 
tion, after which the fluctuating fields are described by a 
random scalar potential y(r, i), rather than by the vector 
potential A(r, t), 



A(r,t)=cJ V r <p(r,t')dt' 



(we assume there are no magnetic fields applied to the 
system). Defining the Cooperon C(r, r'; t, t')) in the new 



gauge by the relation 



C(r,r';t > f)=C(r,r';t,t')) 

( ft/2 ft'/2 

x exp I ie / ip(r, t")dt" - ie <p(r', t")dt' 



(17) 



-ie 



-t/2 1'/2 

<p(r,t")dt" -ie / v {r\t")dt'' 



we obtain the equation 



d 



d 2 



- + '~<p(r, t/2) - -<p{r, -t/2) - D— C(r, r'; t, t') = S (r - r>) 8{t f ) 



at 



dr 2 



(18) 



r 



Note that C(r,r;t,—t) = C(r,r;t,—t), and thus C can 
be used instead of C in for evaluation of the WL correc- 
tion !tT5|) . 

Returning to the consideration of a granular array, 
we assume that the intra-grain conductance is high, 
.9gr ^ 5t- Then the fluctuating potential tp(r,t) does 
not vary from point to point within a single grain, while 
exhibiting random fluctuations of the inter-grain poten- 
tial differences. This allows us to coarse-grain function 
(p(r, i), replacing its dependence on r by the dependence 
on the grain number n. We also can simplify the spatial 
dependence of the Cooperon C(r, r'; t, £')), in case we are 
interested in times long compared to the intra-grain dif- 
fusion time. Indeed, in that case C does not change while 
r or r' vary within a grain. Therefore, the dependence 
of the Cooperon on the coordinates may be replaced by 
the dependence on the grain numbers n and n' which the 
coordinates r and r' belong to. The resulting coarse- 
grained equation for the Cooperon reads: 



d ie , , . ie , , . otS , 
- + - v>„(*/2) - -^Vn{-t/2) + N^- \ C nn ,(t,t') 



g T 6 



— Y J C knl {t,t') = 5 nnl 8{t-t'). 

IT • ■ 



(19) 



Here N is the number of junctions to a single grain (i.e., 
the coordination number of the lattice of grains; N = 4 



for a two-dimensional square lattice), and the summation 
in the last term on the left-hand side runs over N nearest 
neighbors k of the grain n. 

Equations (|15fl and (|19|) provide a convenient starting 
point for evaluation of the weak localization correction at 
temperatures T« T', see Eq. (@J). At higher tempera- 
tures, the spatial dispersion of the fluctuating potentials 
and of the Cooperon inside a grain becomes important. 

The temperature domain T < T* is separated in 
two characteristic regions by the scale T cr , Eq. ©. At 
T <C T" C r, the dependence of Cooperon C on n — n' is 
smooth, and the finite difference equation (|19f) can be re- 
placed by the corresponding differential equation, which 
essentially returns one to the continuous-medium case, 
see Eq. (|18H . Weak localization corrections in this case 
are studied in detail in Refs. IMla. Below we concentrate 
on the temperatures above the crossover. 



IV. QUANTUM CORRECTION TO 
CONDUCTIVITY ABOVE THE CROSSOVER 
TEMPERATURE 



In the temperature regime of interest, 

T* > T > T cr , (20) 
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as we have explained in the Introduction, electron tra- 
jectories are classified according to the number of tunnel 
junctions they cross — the longer are the trajectories, the 
less significant are their contributions. It means that the 
matrix C nn i rapidly decays away from the diagonal. The 
biggest matrix elements are C nn , and the most important 
trajectories are those which do not leave the grain. Eq. 
(|15fl implies that these trajectories do not contribute to 
the weak localization correction, and one needs to con- 
sider the next-order contribution coming from trajecto- 
ries crossing a single junction once. This leads us to 
Eq. © and also allows us to verify the existence of the 
crossover temperature Eq. J2J). 

At T ^ T cr we expect strong fluctuations of the po- 
tential differences between the grains, making coherent 
returns of an electron to the grain of its origin improb- 
able. The returns are described by the term in Eq. i|19|) 
containing the sum over k. Neglecting that term, we find 
for the diagonal component C nn {tt') of the Cooperon 

q£J?(i. = *(* - ex p [~ N ^( t - 

ft/2 r-t/2 ) 

-ie (p n (t")dt" -ie y n {t")dt" \ . (21) 

Jt'/I J-t'/2 j 

The phase factors here reflect the specific gauge we used 
in Eq. (JTOJ). 

Next, we write the Cooperon equation i|19fl for the 
nearest-neighbor sites m and n, 



ie ie S * 

— + -(pmit/2) - ~(p m (-t/2) + JV^L J- C mn {t, t') 



tation, we obtain 



(22) 



The terms with C n i >n describing the grains n' separated 
from n by two tunnel junctions, are small and can be 
omitted in this approximation. Using Eqs. 1)21(1 and l)22[l. 
we obtain for the neighboring grains m and n 

C mn (t, = 9 -^-B{t - t') exp \~N 9 -^{t - t')\ £ dh 

{rti/2 r-ti/2 
—ie / (f n (t")dt" — ie / ip n (t")dt" 
Jt'/2 J-t>/2 

rti/2 p- ii/2 "I 

+ie / tp m (t")dt" + ie ip m (t")dt" \ . (23) 

Jt/2 J-t/2 j 

This expression has to be averaged over the Gaussian 
fluctuations of the field <p. Its correlation function is de- 
fined by the fluctuation-dissipation theorem and reads^ 



''(ipip)(q,uj) = -Im 



2T 1 
uj n(q,w)' 



where II is the polarization operator, calculated for a 
granular medium in Ref. Ildi In the space-time represen- 



e 2 {<p n (t)<p n ,{t')) 



TtTd 2 



91 



5(t - 1') 



dq 



iq(n-n') 



(2^) d Ea( 1 -COsM)' 



(24) 



where the summation in the denominator is carried over 
all available Cartesian components a = x,y, and q a are 
the basis vectors of square lattice of grains. 

Performing the averaging in Eq. with the help of 
Eq. Q24JI is cumbersome but straightforward, since for 
Gaussian fields (exp(z . . . )) = exp(— (. . . )/2). For the 
weak localization correction, we obtain 



(5(7WL 



-A 



9tS 
T ' 



(25) 



.4 



dq 



1 — cos q a d 



N 2 V^\J (2^) d ^ Q (l-cos ga d) 



where V is the volume of the grain (V = d 2 in 2D). Note 
that the interference correction Eq. I|25|l depends on tem- 
perature and on the type of lattice the grains form. The 
dependence on the lattice type comes through the coef- 
ficient A; for a square 2D lattice we find A = 1/4. One 
can easily generalize the evaluation of A onto the case 
of a triangular and more complicated lattices, eventually 
even describing disordered media like ceramics. 



V. MAGNETIC FIELD EFFECT 

As the estimate Eq. suggests, the action of the 
magnetic field on Cooperon is two-fold. A part of the 
Cooperon suppression comes from the intra-grain elec- 
tron motion, and another part stems from the magnetic 
field effect on the inter-grain coherence. Since g gr S> gT, 
the interesting range of the fields corresponds to a small 
flux penetrating a grain, Bd 2 <C $o- We can then con- 
sider the effect of magnetic field on the Cooperon within 
a grain perturbatively. To implement the perturbation 
theory, it is convenient to use a "tailored" to the grains 
shape gauge of the magnetic field B. For dcfinitcncss, 
we concentrate on the case of a two-dimensional array of 
"flat" grains connected by point-like tunneling contacts, 
see Fig. 1. We define the gauge for the points within the 
grains by the relations 



A (r) = rx V^(r)+Ai 



B, ip„(r e b n ) = 0. 

(26) 

Here r is the normal to the plane of the grains, and 
b n is the boundary of the n-th grain. The second and 
third relations in Eq. (|26|l fully define the boundary prob- 
lem for a scalar function ip n (r). The constants A n are 
tuned in such way that the vector-potential is continu- 
ous at the points of contact between the grains. Up to 
a discrete analogue of the gradient of a scalar function, 
these constants are determined fully by the solution of 



7 



the boundary problems for all tfj m (r). It is clear that the 
discrete version of curl applied to A n must be equal B 
upon averaging over the array; the characteristic differ- 
ence A n — A m for two nearby junctions is of the order 
A n A m ^ B • d. 

In the definition of the Cooperon, it is convenient to 
present again the coordinates as pairs {r, n} and {r', n'} 
which point explicitly to the label of grains the two points 



r,r' belong to. In addition, we multiply the Cooperon 
defined in Eq. I|17|l by yet one more gauge factor, 

C mn (r,r';t,t') = C^ n (r,r';t,t') exp(iA m ■ r - iA n ■ r'), 

. (27) 

In these new notations, the equation for Cooperon in the 
absence of tunneling has the form 



^ + ^(r, t/2) - ^(r, -t/2) - D gl . (J- - i-r x VW»(r)) } Ct n {r, r>; t, t>) = S mn S (r - r') 6(t - t'), (28) 

I 

where D gI ~ d 2 g gr S is the diffusion coefficient within the with the unit quantum, Bd 2 -C w ^ may treat the 

grain (here g gr 5 is the Thouless energy for the electron effect of a magnetic field within a grain perturbatively. 

motion within a grain). With the defined gauge Eq. (|2f)|) . Considering the low-energy limit, T <C g gr S, and taking 

the normal to the boundary component of A(r) is zero. into account the boundary conditions for C^, we start 

Thus, the magnetic field does not affect the boundary perturbations from r-independent Cooperon C^ n (f, t'). 

conditions for Cooperon, i.e. the normal component of In the presence of inter-grain tunneling, the correspond- 

dC^/dr at the boundary is zero. ing generalization of Eq. I|I9|I reads 
As long as the flux piercing one grain is small compared 




(Bd 2 



~M*/2) - ^Vm(-*/2) + N^- Ct n (t,t') 



£ e^M-^)ct n (t, t>) = S mn S(t - t'). 

k 

I 



(29) 



Here the magnetic field dependence 



ag gI 6 



Bd 2 



D e 2 
d 2 c 2 



d 2 r\VMr)\' 



comes from the V'n-dependent term in Eq. I|28J) integrated 
over the volume of a single grain; a ~ I is the dimen- 
sionless coefficient depending on the grains shapes. The 
vector rkn points to the junction between grains k and 



The discreteness of the medium is adequately ac- 
counted for by the structure of Eq. jSSjl. However, 
the discreteness is not important in the domain of low 
temperatures and relatively low fields, T -C T cr and 
B <C i?0 S . There we can replace the left-hand side of 
Eq. H29fl by its gradient expansion. After the expansion 
and replacement of the grain number n by the corre- 



sponding coarse-grained coordinate R, we find 



d [Bd 2 



D 



V 



R 



-A(R) 



(30) 



+^(R, t/2) ~ l ^(R, -t/2)} C*(R, R'\ t, t') 
= S(R-R')S(t-t'). 

The second term in the left-hand side here reflects the 
suppression of interference by the magnetic flux pene- 
trating the grains. Apart from that term and from the 
value of the effective diffusion constant D = ■n^ 1 g^5d 2 , 
which reflects the granularity of the medium, this equa- 
tion is identical to that of a homogeneous thin film. Using 
the known results^ for the films, we find the magnetore- 
sistance of a granular array, 



5(Tmr{B,T) = 5(j W l{B,T) - Sa W L(0,T) 



(31) 
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FIG. 3: Magnetic field dependence (1311 of the conductivity at 
low temperatures, T <C T^gT / <?cr. The horizontal axis here 
is proportional to the logarithm of the applied magnetic field, 
y = ln(_Bd 2 /$o) + ln(g gI /gT); the vertical axis is the magne- 
toresistance in dimensionless units, <5<7mr = (e 2 /2h)f(y) with 
f(v) — V + m (l + e!/ )i see Eq. 1311 . The crossover between 
ln_B and 21nB dependences is clearly seen. 

(we dispensed with the factor a ~ 1 here) . The two field 
scales introduced in Eqs. JSJ and Q can be obtained 
from a comparison [in the argument of logarithmic func- 
tion Eq. JSU] °f the dephasing term T/T cr with the linear 
and quadratic in B terms, respectively. At lowest tem- 
peratures, there is a clear crossover in the 5<jwl vs. B 
dependence from Sayjh oc InB to <5<7wl cx 21ni3. Note 
that the crossover occurs in the 2D regime, where the 
typical closed path for a coherent electron motion spans 
many grains. It is remarkable that even in the 2D regime 
there is a clear difference in the magnetoresistance of a 
granular system from that of a homogeneous film, see 
Fig. 3. 

VI. DISCUSSION 

Let us now discuss the temperature dependence of 
the conductance in various magnetic fields (Fig. 4). In 
zero field, the weak localization correction behaves as 
<5cwl/co ~ g^ bxT at T < T CI and then crosses over 
to the power-law behavior, £<7wl/0u ~ T a ./(jg^T'), at 
higher temperatures. Finite magnetic field leads to sup- 
pression of the WL correction even at the lowest tem- 
perature. Thus, at B < ($ /d 2 )(g T /gg r ) the WL cor- 
rection becomes temperature- independent. (Note that 
($o/ d 2 ){gT I ' 9%v) <C B s ^ .) At higher temperatures, the 
dimensionless conductivity <5ctwl/ f"o has the same tem- 
perature dependence as at B = 0. In higher fields, 
(4>o/<i 2 )((?T/9gr) <C B <C B s ^, the same low-temperature 
saturation occurs at T = (Bd 2 /$o) 2 (g gT /gT)T cr , see 
Fig. 4. In the highest fields, B > B s ?, the WL correc- 
tion is suppressed for all trajectories — even those lying 
within a single grain, and the WL correction disappears 
at all temperatures. 

Note that all magnetic fields which we have discussed 
above are too small to change orbital dynamics of elec- 
trons. Indeed, the cyclotron radius, r c = mvFc/eB must 
be smaller than the mean free path I in order to affect 



G 




Tor T 

FIG. 4: Sketch of the temperature dependence of the con- 
ductance in various magnetic fields: B — (1); B <C B^ 
(2). The temperatures at which curve 2 departs considerably 
from curve 1 depend on the applied field; these temperatures 
are of the order of T cr (Bd 2 /$ ) and T CT {Bd 2 /$ ) 2 (s g r / 'g T ) 
for B < ($ /d 2 )(W£fT) and B > (®o/d 2 )(g gI /g T ), respec- 
tively. 



the electron motion. This corresponds to magnetic fields 
B > ($o/d 2 )(T£) _1 , with t being the momentum re- 
laxation time in a grain. Since t5 <C 1 (conditions for 
metallic diffusive behavior), such fields are well outside 
our consideration range. 

Apart from the weak localization correction, there is 
one more temperature dependent contribution to the con- 
ductance — interaction correction. For granular media, 
it was calculated for all temperatures in Ref. El It 
crosses over from low- to high-temperature regime at the 
temperature g^S, which is different (much lower) than 
T cr . For a two-dimensional array, this correction is log- 
arithmic at any temperatures; for T 3> <7t<5, one has 
Sa/ao ~ g^}\Yi(g^Ec /T), where Eq is the charging en- 
ergy in a single grain. The temperature dependence of 
the interaction correction is featureless at T ~ T cr , and 
therefore it should not mask the crossover in the temper- 
ature dependence of the WL correction. The interaction 
correction is also independent of the magnetic field. Thus 
it does not affect the crossover in the magnetic field de- 
pendence of the conductance, which is induced by the 
granular structure. The measurements of the conduc- 
tance therefore can be used to characterize the medium. 

Let us finally give some estimates. We consider metal- 
lic grains of a size of 500 nm, which can be easily pro- 
duced lythographycally-. They have the level spacing 
of order ~ 20mK. Choosing <?t = 10, we obtain 

the crossover temperature T CT = 2K, that can be easily 
observed experimentally. 
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Note added: After completing this work, we noticed 
that a formula similar to our Eq. was derived, with 



a different method, in the work by Biagini et al very 
recently^. We are grateful to Andrei Varlamov for the 
discussion of relation between the two works. 
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